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1. Introduction 

To understand the origin of high-energy eosmie rays is one of the open key questions in as- 
tropartiele physies [1,2]. An inspiring artiele, published in 2003 [3] heralded the renaissanee of radio 
deteetion of extensive air showers with the ultimate goal to measure the properties of eosmie rays 
with this teehnique and the pioneer experiments LOPES [4] and CODALEMA [5] were initiated. 
The big sueeess of these pathfinders stimulated further investigations of the radio emission of air 
showers on larger seales, with installations sueh as Tunka-Rex [6], AERA (Auger Engineering Ra¬ 
dio Array) at the Pierre Auger Observatory, and the EOEAR radio teleseope. Signifieant progress has 
been aehieved in the last deeade [7] and we now understand the emission proeesses of the radio waves 
in the atmosphere. Most of the emission is due to the interaetion of the shower with the magnetie field 
of fhe Earfh, whieh leads fo a fransverse eurrenf in fhe shower. In addifion fo fhis emission, fhe over- 
abundanee of eleefrons in fhe shower fhaf are eolleefed from afmospherie moleeules leads fo a eurrenf 
in fhe direefion of fhe shower. 

2. Radio Air-Shower Detectors 

Two modern radio defeefors for exfensive air showers are fhe EOEAR radio feleseope (in parfie- 
ular, ifs dense eore in fhe Nefherlands) and AERA af fhe Pierre Auger Observafory in Argentina. The 
layouf of fhe fwo sef-ups is drawn fo seale in Pig. 1. 

LOFAR radio telescope The LOPAR key seienee projeef Cosmie Rays represenfs one of fhe six 
seienfifie key objeefives of EOEAR [8,9]. EOEAR is a digifal radio feleseope. Ifs antennas are spread 
over several European eounfries and are used fogefher for inferferomelrie radio observations in fhe 
frequeney range of 10 - 240 MHz. The density of anfennas inereases towards fhe eenfer of EOEAR, 
whieh is loeafed in fhe Nefherlands. Here, abouf 2400 anfennas are elusfered on an area of roughly 
10 km^. This high densify of anfennas makes EOEAR fhe perfeef fool fo sfudy fealures of fhe radio 
emission ereafed by exfensive air showers. Air-shower measuremenfs are conduefed based on a frigger 
reeeived from an array of seinfillafors (EORA) [10,11], whieh resulfs in a read-ouf of fhe ring buffers 
fhaf sfore fhe raw volfage fraees per antenna for up fo 5 s. EOEAR eomprises fwo fypes of anfennas, 
reeording radio emission in low-frequeney band from 10 fo 90 MHz and also in fhe high-frequeney 
band (110-240 MHz). 
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Fig. 1. Layout of AERA at the Pierre Auger Observatory and the dense core of LOFAR - drawn to scale. 


AERA at the Pierre Auger Observatory AERA is a radio extension of the Pierre Auger Observa¬ 
tory. The Observatory is a 3000-km^ hybrid cosmic-ray air-shower detector in Argentina [12,13], with 
an array of 1660 water-Cerenkov detectors and 27 fluorescence telescopes at four locations on the pe¬ 
riphery [14]. The area near the Coihueco fluorescence detector contains a number of low-energy en¬ 
hancements, including AERA. AERA is located in a region with a higher density of water-Cerenkov 
detectors (on a 750 m grid) and within the field of view of HEAT [15], allowing for the calibration 
of the radio signal using super-hybrid air-shower measurements, i.e., recording simultaneously the 
fluorescence light, the particles at the ground, and the radio emission from extensive air showers. 

Since March 2015 AERA consists of 153 autonomous radio-detection stations, distributed with 
different spacings, ranging from 150 m in the dense core up to 750 m, covering an area of about 
17 km^. Different types of antennas are used, including logarithmic periodic dipoles and active bowtie 
antennas (named ’’Butterfly antennas” [16]) covering the frequency range from 30 to 80 MHz [17,18]. 

3. Precision measurement of the radio emission in air showers 

EOEAR combines a high antenna density and a fast sampling of the measured voltage traces in 
each antenna. This yields very detailed information for each measured air shower and the properties 
of the radio emission have been measured with high precision. At the Pierre Auger Observatory air 
showers are measured simultaneously with various detector systems: radio detectors, fluorescence 
light telescopes, water-Cerenkov detectors, and underground muon detectors. This unique combi¬ 
nation yields complementary information about the showers and allows to investigate correlations 
between the various shower components. Some important aspects of radio emission in air showers 
are reviewed in the following. We focus on radio emission in the frequency range 30 - 80 MHz, only 
one result (Pig. 3 right) deals with higher frequencies. 
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Fig. 2. Left: Footprint of an air shower measured with AERA at the Pierre Auger Observatory, the boxes 
representing the radio-detector stations. The colored background represents predictions of the radio signal ac¬ 
cording to simulations [19]. Right: The arrival time of the signals, relative to aplane as a function to the distance 
to the shower axis as measured with LOFAR. The lower right graph illustrates the arrival time differences with 
respect to a fit (hyperboloid) [20]. 


Lateral distribution function of the radio signals The footprint of the radio emission recorded at 
ground level is not rotationally symmetric [19,21,22], such as, e.g., the particle content of a shower, 
see Fig. 2 (left). Radio emission is generated through interactions with the Earth magnetic field, which 
yields a bean-shaped footprint on the ground. The correct reference system is in the shower plane (per¬ 
pendicular to the shower axis), with one coordinate axis along the vxS direction and the second axis 
along the v x (v x d) direction. Here v is the propagation velocity vector of the shower (parallel to the 
shower axis) and B represents the direction and strength of the Earth magnetic field. The measured 
power in the frequency range 30 - 80 MHz is plotted as a function of the distance to the shower 
axis in Eig. 3 (left). Eor example at a distance of 200 m from the shower axis, ambiguities are visible 
in this one-dimensional projection: the recorded signal strength is a function of the azimuth angle, 
which results in the visible structure. Historically, the lateral distribution of the radio signal on the 
ground has often been parameterized with a simple exponential function (e.g., [23]). However, the 
EOEAR measurements suggest that the radio emission should be parameterized by a more complex 
expression: a two-dimensional Gaussian function is used to describe the approximately exponential 
fall-off at large distances from the shower axis. A second (smaller) two-dimensional Gaussian func¬ 
tion is subtracted from the first one to describe the ring structure of the signal close to the shower 
axis. To reproduce the observed bean shape, the centers of both Gaussian functions are slightly offset. 
The power at position ix',y') in the shower plane (perpendicular to the shower axis) is described as 
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Where the parameters are the scaling factors A+ and A_ (where in general A+ >> A_), the width 
of the Gaussian functions a-+ and cr_, and the centers of the two-dimensional Gaussian distributions 
(X+, F+) and (A_, F_), respectively. 

First quantitative measurements in the frequency range 120-240 MHz Radio emission from 
extensive air showers has also been recorded with the EOEAR high-band antennas in the 200 MHz 
frequency domain [24]. The measured power is depicted in Eig. 3 (right) as a function of the distance 
to the shower axis. A clear maximum is visible at distances around 120 m in this one-dimensional 
projection, indicating a clear (Cerenkov) ring structure. Such rings are predicted from theory: rela- 

















0.40 



-50 0 50 100 150 200 250 300 350 400 0 100 200 300 


Distance to shower axis [m] distance (m) 

Fig. 3. Total power measured in an air shower as a function of the distance to the shower axis (in the shower 
plane) as measured by LOFAR in the frequency range 30 - 80 MHz (left) [22] and 110-190 MHz (right) [24]. 


tivistic time compression effects lead to a ring of amplified emission, which starts to dominate the 
emission pattern for frequencies above ~ 100 MHz. The LOFAR data clearly confirm the importance 
to include the index of refraction of air as a function of height into calculations of the radio emission. 
Shape of the shower front The precise shape of the radio wavefront is a long-standing issue. 
In the literature different scenarios have been discussed: a spherical, conical, or hyperbolic shape 
(e.g. [25]). The LOFAR findings clearly indicate that a hyperboloid is the best way to describe the 
shape of the measured wavefront [20]. A hyperboloid asymptotically reaches a conical shape at large 
distances from the shower axis and can be approximated as a sphere close to the shower axis. A 
measured wavefront of a shower registered with LOFAR is shown in Fig. 2 (right). The time difference 
relative to a plane is plotted as a function of distance to the shower axis. The line indicates a fit of 
a hyperboloid to the measured data. The lower part of the graph shows the time differences of the 
individual antennas with respect to the fit function. 

Polarization of the radio signal The radio emission in extensive air showers originates from dif¬ 
ferent processes. The dominant mechanism is of geomagnetic origin [4,26]. Electrons and positrons 
in the shower are accelerated in opposite directions by the Lorentz force exerted by the magnetic field 
of the Earth. The generated radio emission is linearly polarized in the direction of the Eorentz force 
(v X B), where v is the propagation velocity vector of the shower (parallel to the shower axis) and B 
represents the direction and strength of the Earth magnetic field. A secondary contribution to the radio 
emission results from the excess of electrons at the front of the shower (Askaryan effect) [27,28]. This 
excess is built up from electrons that are knocked out of atmospheric molecules by interactions with 
shower particles and by a net depletion of positrons due to annihilation. This charge excess contribu¬ 
tion is radially polarized, pointing towards the shower axis. The resulting emission measured at the 
ground is the sum of both components. Interference between these components may be constructive 
or destructive, depending on the position of the observer/antenna relative to the shower. The emission 
is strongly beamed in the forward direction due to the relativistic velocities of the particles. Addition¬ 
ally, the emission propagates through the atmosphere, which has a non-unity index of refraction that 
changes with height. This gives rise to relativistic time-compression effects, most prominently result¬ 
ing in a ring of amplified emission around the Cerenkov angle, see Eig. 3. By precisely measuring 
the polarization direction of the electric field at each antenna position, the ratio a between the con¬ 
tributions from the Askaryan effect and the geomagnetic emission is measured. The resulting ratios 
as measured at the Auger Observatory for individual air showers are shown in Eig. 4 (left) [29]. An 
average value of a = 14% + 2% is obtained. At EOEAR the ratio a has been measured as a function of 
the distance to the shower axis for showers with different zenith angles, as depicted in Eig. 4 (right) 
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Fig. 4. Radio emission processes in the atmosphere. The ratio a between contributions due to the Askaryan 
effect and the geomagnetic emission is plotted as measured at the Auger Observatory for different air showers 
(left) [29] and as a function of the distance to the shower axis (in the shower plane) for showers with different 
zenith angles as measured by LOFAR (right) [30]. 


[30]. The figure illustrates that the contribution through the Askaryan effect increases with increasing 
distance to the shower axis and it is more pronounced for vertical showers (with small zenith angle). 
Confirmation of simulation codes The detailed investigations of the properties of the radio emis¬ 
sion and the comparison between measurements and predictions from simulations demonstrate that 
the simulation code CoREAS [31] fully describes all relevant features of the radio emission in air 
showers and the predictions can be used to interpret the measured air-shower data. 

Probing atmospheric electric fields during thunderstorms Radio detection of air showers can 
also be used for auxiliary science, such as the measurements of electric fields in the atmosphere during 
thunderstorms [32-34]. The footprint of the radio emission from an air shower, which developed 
during a thunderstorm is shown in Fig. 6 (right) [34]. The intensity and polarization patterns of such 
air showers are radically different from those measured during fair-weather conditions. The figure 
illustrates the polarization as measured with individual LOFAR antennas (arrows) in the shower plane. 
An arrow labeled ’’normal” indicates the expected polarization direction for fair-weather conditions. 
LOFAR antennas are grouped into circular stations, of which seven are depicted. The position of 
the shower axis, orthogonal to the shower plane, is indicated by the intersection of the dashed lines. 
With the use of a simple two-layer model for the atmospheric electric field, these patterns can be 
well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study 
atmospheric electric fields. 

4. Measuring properties of cosmic rays with the radio technique 

The ultimate objective is to fully characterize the incoming cosmic ray with the radio measure¬ 
ments, i.e., to determine its arrival direction, its energy, and the particle type/mass (see, e.g., [35]). 
Arrival direction A precise description of the shape of the shower front is essential to correctly 
reconstruct the direction of the incoming cosmic ray. Different shapes have been investigated at LO- 
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Fig. 5. Determining the properties of cosmic rays with radio measurements [22]. Left: Energy of the cosmic 
ray (determined from simulations and measured with the LORA air shower array) as a function of the integrated 
measured radio power. Right: Width of the footprint of the radio emission as a function of the distance to the 
shower maximum. 


FAR as discussed above [20]. Reconstructing the arrival directions of the same measured showers 
with different wavefront shapes results in differences in fhe reconsfrucfed arrival directions of typi¬ 
cally less than 1° between a plane and a hyperboloid and typically less than 0.1° between a cone and 
a hyperboloid. Based on these investigations it is expected that the angular resolution for the arrival 
direction of the shower is better than 1°. Corresponding analyses at the Auger Observatory are under 
way, expecting a similar precision. 

Energy The parameters of the function to model the intensity pattern of the radiation on the ground 
(as described above) are sensitive to the properties of the shower-inducing cosmic rays [22]. The 
integral of the measured power density is proportional to the shower energy. This is illustrated in 
Fig. 5 (left). The shower energy is plotted as a function of a parameter, that is proportional to the 
integrated power. The measured signal strength is corrected for a factor, which depends on the angle 
between the shower axis and the direction of the geomagnetic field. This is necessary, since fhe 
dominant geomagnetic emission strongly depends on this angle. The shower energy is determined 
in two ways in the figure: if is derived from Monfe Carlo simulafions (using CORSIKA [36] and 
CoREAS) and if has been measured wifh fhe particle defector array af LOFAR (LORA). A clear 
correlafion can be seen befween fhe shower energy and fhe measured radio infensify on fhe ground. 

A very sophisficafed analysis is ongoing af fhe Auger Observatory [37-39]. The energy content of 
the radio emission in air showers in the frequency range 30 - 80 MHz has been measured on an abso¬ 
lute scale. A 1 EeV air shower, arriving perpendicularly to the geomagnetic field radiates 15.8 MeV in 
the radio frequency regime. Since radio emission is not suffering from affenuafion in fhe afmosphere 
and fhe radiafion sfrengfh can be calculafed from firsf principles, measuring fhe radio emission seems 
fo be a very promising fool fo esfablish fhe absolufe energy scale of air-shower defecfors. 

Depth of the shower maximum Several methods are studied at LOEAR and the Auger Observatory 
to derive the type of the incoming cosmic ray or its atomic mass A from the radio measurements. The 
key is to measure the depth of the shower maximum in the atmosphere, which is proportional to In A. 
The width of the measured footprint is proportional to the geometrical distance from the antennas to 
the position of the shower maximum. This correlation is shown in Eig. 5 (right). This method is used 
at LOEAR and at the Auger Observatory to estimate the cosmic-ray mass [19,22]. 

To obtain a precise value of the depth of the shower maximum for each shower, we apply the 
following procedure, originally developed at LOEAR and now also applied at the Auger Observatory 
[19,40]. The shower direction and energy are obtained from the measured signals in the particle de- 
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Fig. 6. Left: Determining the depth of the shower maximum from the radio measurements. The;^^^ values of 
simulated showers are shown as a function of the depth of the shower maximum [40]. Right: Footprint of the 
radio emission of an air shower, which developed during a thunderstorm [34]. 


lectors and the radio antennas. With these parameters simulations are initiated, simulating the devel¬ 
opment of air showers and the accompanying radio emission, using the CORSIKA and the CoREAS 
codes, for primary protons and iron nuclei. Due to fluctuations in the individual particle interactions 
in an air shower this results in cascades with a wide distribution for the depth of the shower maximum 
(all of them having the same energy and direction of incidence as the measured one). The predicted 
signals in the particle detectors and the radio antennas are then compared to the measurements. A 
method is used to determine the simulated shower which best matches the measured values. This 
yields a value for the depth of the shower maximum. The method is illustrated in Fig. 6 (left). The 
values are depicted as a function of the depth of the shower maximum X^ax- A parabola is fitted 
to the x^ values and the minimum of this function gives the estimated Xmajcvalue for the measured 
shower. With this method X-axis determined with an accuracy of 17 g/cm^ at LOFAR. At the Auger 
Observatory this method is used to compare the depth of the shower maximum as obtained from 
the radio measurements with the X^oxValues measured with the fluorescence detectors. Quantitative 
analyses of the achieved accuracy are under way. 

At the Auger Observatory other mass-sensitive observables are also investigated: the shape of 
the shower front is sensitive to the geometrical distance to the shower maximum; the shape of the 
shower front is approximated with a hyperboloid and the angle between the shower plane and the 
hyperboloid (at large distances from the shower axis) is sensitive to X^ax- The spectral shape of the 
measured radio signal is also sensitive to X^ax [41]- This could be a very promising method for a 
sparse radio array, where the only measurable radio signals are detected in a small number of radio 
detectors. 
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